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Abstract. We prove the following estimate for the spectrum of the normalized 
Laplace operator A on a finite graph G, 

1 - (1 - k[t])* < Xi < ■ ■ ■ < X N -i < 1 + (1 - k[t])*, V integers t > 1. 

Here k[t] is a lower bound for the Ollivier-Ricci curvature on the neighborhood graph 
G[t] (here we use the convention G[l] — G), which was introduced by Bauer- Jost. 
In particular, when t = 1 this is Ollivier's estimate k < Ai and a new sharp upper 
bound \n-i < 2 — fc for the largest eigenvalue. Furthermore, we prove that for any 
G when t is sufficiently large, 1 > (1 — k[t])i which shows that our estimates for Ai 
and Ajv_i are always nontrivial and the lower estimate for Ai improves Ollivier's 
estimate k < Ai for all graphs with k < 0. By definition neighborhood graphs 
possess many loops. To understand the Ollivier-Ricci curvature on neighborhood 
graphs, we generalize a sharp estimate of the curvature given by Jost-Liu to graphs 
which may have loops and relate it to the relative local frequency of triangles and 
loops. 



1. Introduction 

In this paper, we utilize techniques inspired by Riemannian geometry and the the- 
ory of stochastic processes in order to control eigenvalues of graphs. In particular, 
we shall quantify the deviation of a (connected, undirected, weighted, finite) graph 
G from being bipartite (a bipartite graph is one without cycles of odd lengths; equiv- 
alently, its vertex set can be split into two classes such that edges can be present 
only between vertices from different classes) in terms of a spectral gap. The operator 
whose spectrum we shall consider here is the normalized graph Laplacian A (it is 
unitarily equivalent to the one studied in Chung jH]). This is the operator underlying 
random walks on graphs, and so, this leads to a natural connection with the the- 
ory of stochastic processes. We observe that on a bipartite graph, a random walker, 
starting at a vertex x at time and at each step hopping to one of the neighbors of 
the vertex where it currently sits, can revisit x only at even times. This connection 
then will be explored via the eigenvalues of A. More precisely, the largest eigenvalue 
\n-i of A is 2 iff G is bipartite and is < 2 else. Therefore, 2 — Xn-i quantifies the 
deviation of G from being bipartite, and we want to understand this aspect in more 
detail. In more general terms, we are asking for a quantitative connection between 
the geometry (of the graph G) and the analysis (of the operator A, or the random 
walk encoded by it). Now, such connections have been explored systematically in 
Riemannian geometry, and many eigenvalue estimates are known there that connect 
the corresponding Laplace operator with the geometry of the underlying space M, 
see e.g. Li-Yau [Hj, Chavel [I]. The crucial role here is played by the Ricci curvature 
of M. In recent years, a kind of axiomatic approach to curvature has been developed. 
This approach encodes the abstract formal properties of curvature and thereby makes 
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the notion extendible to spaces more general than Riemannian manifolds. By now, 
there exist many notions of generalized curvature, and several of them have found 
important applications, see Sturm [201121], Lott-Villani [16], Ollivier [IS] . Ohta [17] . 
Bonciocat-Sturm [2], Joulin- Ollivier [12] and the references therein. The curvature 
notion that turns out to be most useful for our purposes is the one introduced by 
Ollivier [IS]. One of our main results says that 



(1.1) X N ^<2-k 

where k is a lower bound for Ollivier's Ricci curvature. This matches quite well with 
Ollivier's result that the smallest nonzero eigenvalue Ai satisfies 

(1.2) Ai > k. 

In fact, one of the main points of the present paper is to relate such upper and lower 
bounds via random walks. In fact, as in Bauer- Jost [I], we translate this relationship 
into the geometric concept of a neighborhood graph. The idea here is that in the 
t-th neighborhood graph G[t] of G, vertices x and y are connected by an edge with a 
weight given by the probability that a random walker starting at x reaches y after t 
steps times the degree of x. We note that even though the original graph may have 
been unweighted, the neighborhood graphs G[t] are necessarily weighted. In addition, 
they will in general possess self-loops, because the random walker starting at x may 
return to x after t steps. Therefore, we need to develop our theory on weighted 
graphs with self-loops even though the original G might have been unweighted and 
without such loops. Since Ollivier's curvature is defined in terms of transportation 
distances (Wasserstein metrics), we can then use our neighborhood graphs in order to 
geometrically control the transportation costs and thereby to estimate the curvature 
of the neighborhood graphs in terms of the curvature of the original graph. As it turns 
out that lower bounds for the smallest eigenvalue of G[2] are related to upper bounds 
for the largest eigenvalue of G, we can then use methods for controlling the smallest 



eigenvalue to also obtain ( 1.1 ). For controlling the smallest eigenvalue, besides Ollivier 



[18] . we also refer to Lin-Yau [15] and Jost-Liu [IT]. In particular, in the last paper, we 
could relate Ai to the local clustering coefficient introduced in Watts-Strogatz [2"3] . 
The local clustering coefficients measures the relative local frequency of triangles, 
that is, cycles of length 3. Since bipartite graphs cannot possess any triangles, this 
then is obviously related to our question about quantifying the deviation of the given 
graph G from being bipartite. In fact, in Jost-Liu [IT], this local clustering has been 
controlled in terms of Ollivier's Ricci curvature. Thus, in the present paper we are 
closing the loop between the geometric properties of a graph G, the spectrum of its 
graph Laplacian, random walks on G, and the generalized curvature of G, drawing 
upon deep ideas and concepts originally developed in Riemannian geometry and the 
theory of stochastic processes. 



2. The normalized Laplace operator, neighborhood graphs, and 

Ollivier-Ricci curvature 

In this paper, G = (V, E) will denote an undirected, weighted, connected, finite 
graph of iV vertices. We do not exclude loops, i.e., we permit the existence of an 
edge between a vertex and itself. V denotes the set of vertices and E denotes the 
set of edges. If two vertices x, y G V are connected by an edge, we say x and y are 
neighbors, in symbols x ~ y. The associated weight function w: V x V — > R satisfies 
w xy = w yx (because the graph is undirected) and we assume w xy > whenever x ~ y 
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and w xy = iff x 7^ y\ thus, we do not permit negative weights. For a vertex x G V, 
its degree d x is defined as d x := ^2 yeV w xy . If w xy = 1 whenever x ~ y, we shall call 

the graph an unweighted one. We will also speak of a locally finite graph G = (V, E), 
which is an undirected, weighted, connected graph and satisfies the property that for 
every x G V, the number of edges connected to x is finite. 

2.1. The normalized graph Laplace operator and its eigenvalues. In this 
subsection, we recall the definition of the normalized graph Laplace operator and 
state some of its basic properties. In particular, we will emphasize the relations 
between eigenvalues of the Laplace operator and random walks on graphs. 

Let C(V) denote the space of all real- valued functions on the set V. We attach 
to each vertex x G V a probability measure m x (-) and then can define the Laplace 
operator. 

Definition 1. The Laplace operator A is defined as 

(2-1) Af(x) = f(y)m x (y) - f(x), V/ G C(V). 

y&V 

The measure m x {-) can also be considered as the distribution of a 1-step random 
walk starting from x. We will choose 

f ^ jf y ^ X 

< 2 ' 2 > = { * • ot L™e, 

in the following. Note that x ~ x is possible when x has a loop. With this family 
of probability measures {m 1 (-)} ) A is just the normalized graph Laplace operator 
studied in Grigoryan [TU] and Bauer- Jost [I] and is unitarily equivalent to the Laplace 
operator studied in Chung jS]. 

Remark 1. On a graph G without loops, we can also consider a lazy random walk. 
A lazy random walk is a random walk that does not move with a certain probability, 
i.e. for some x we might have m x (x) 7^ 0. In this case, the lazy random walk on G 
is equivalent to the usual random walk on the graph G Uzy that is obtained from G by 
adding for every vertex x a loop with a weight d x m x (x). 

We also have a natural measure [i on the whole set V, 

H{x) := d x , 

which gives an inner product structure on C(V). 

Definition 2. The inner product of two functions f,g G C(V) is defined as 
(2.3) (f,g) fl = J2f(^9(x)fi(x). 

Then C(V) becomes a Hilbert space, and we can write C(V) = l 2 (V,fi). By the 
definition of the degree and the symmetry of the weight function, we can check that 

• fx is invariant w.r.t. {m x (-)}, i.e. J2 XIE v m x{y)n{x) = fJ>(y), G V; 

• fi is reversible w.r.t. {m x (-)}, i.e. m x (y)fi(x) = m y (x)fi(y), Vx,y G V. 

These two facts imply immediately that the operator A is nonpositive and self-adjoint 
on the space l 2 (V, ji). We call A an eigenvalue of A if there exists some / ^ such 
that 



(2.4) 



A/ = -Xf. 
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Using this convention it follows from the observation that A is self-adjoint and non- 
positive that all its eigenvalues are real and nonnegative. In fact, it's well known that 
(see e.g. Chung [8]) 

(2.5) = A < Ai < • • • < Ajv_i < 2. 

Since our graph is connected we actually have < Ai. In Chung [B] it is shown, by 
proving a discrete version of the Cheeger inequality, that Ai is a measure for how 
easy/difficult it is to cut the graph into two large pieces. Furthermore, Atv_i = 2 if 
and only if G is bipartite. In Bauer- Jost pQ it is shown, by establishing a dual version 
of the Cheeger inequality for the largest eigenvalue, that Xn-i is a measure of the size 
of a locally bipartite subgraph that is connected by only a few edges to the rest in 
the graph. In the following, we will call Ai the first eigenvalue and Aat_i the largest 
eigenvalue of the operator A. 

2.2. Neighborhood graphs. In this subsection, we recall the neighborhood graph 
method developed by Bauer- Jost [lj. 

As discussed above, the Laplace operator underlies random walks on graphs. In 
this section, we discuss the deep relationship between eigenvalues estimates for the 
Laplace operator A and random walks on the graph G by using neighborhood graphs. 

We first introduce the following notation. For a probability measure /i, we denote 



x)m a 



Let 5 X be the Dirac measure at x, then we can write S x P l (-) := S X P(-) = m x (-). 
Therefore the distribution of a t-step random walk starting from x with a transition 
probability m x is 

(2.6) 5 x P\-) = ^2 m x{xi)m xi {x 2 ) ■• •m Xt _ 1 (-) 

x\,...,x t -\ 



for t > 1. In particular, using the measure (2.2), the probability that the random 



walk starting at x moves to y in t steps is given by 

xi,...,x t -i d x d x , d x , , ' 



w *y if / 

4= ' 11 L 



The idea is now to define a family of graphs G[t], t > 1 that encodes the transition 
probabilities of the t-step random walks on the graph G. 

Definition 3. The neighborhood graph G[t] = (V,E[t\) of the graph G = {V,E) of 
order t > 1 has the same vertex set as G and the weights of the edges of G[t] are 
defined by the transitions probabilities of the t-step random walk, 

(2.7) w xy [t] := 5 x P\y)d x . 

In particular, G — G[l] and x ~ y in G[t] if and only if there exists a path of length 
t between x and y in G. 

Remark 2. We note here that the neighborhood graph method is related to the discrete 
heat kernel p t (x , y) on graphs. For more details about the discrete heat kernel see for 
instance Grigoryan [10]. We have 

W X y[t] 



Pt(x,y) 



d x dy 
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Example 1. We consider the following two examples. 

G G[2] 




1/2 

Note that the neighborhood graph G[2] is disconnected. In fact the next lemma shows 
that this is the case because G is bipartite. Note furthermore that E(G) ^ E(G[2\). 

H H[2] 




1 



For this example we have E(H) C E(H[2]). 

These examples shows that the neighborhood graph G[t] is in general a weighted 
graph with loops, even if the original graph G is an unweighted, simple graph. 

Lemma 1. The neighborhood graph G[t] has the following properties (see Bauer- Jost 

W : 

(i) If t is even, then G[t] is connected if and only if G is not bipartite. Further- 
more, if t is even, G[t] cannot be bipartite. 

(ii) Ift is odd, then G[t] is always connected and G[t] is bipartite iff G is bipartite. 

(iii) d x [t] = d x for all x G V . 

Note that {iii) implies that l 2 (V,fi) = l 2 (V,fj,[i\) for all t. In particular, we have 
the same inner product for all neighborhood graphs G[t}. The crucial observation is 
the next theorem. 

Theorem 1 (Bauer- Jost). The Laplace operator A on G and the Laplace operator 
A[t] on G[t] are related to each other by the following identity: 

id - (id- A)' = A[t]. 

The nonlinear relation between A and A [t] leads directly to the following eigenvalue 
estimates. 

Theorem 2 (Bauer- Jost). Let A[t] be a lower bound for the eigenvalue Xi[t] of A[t], 
i.e., Ai[i] > A[t}. Then 

(2.8) 1 - (1 - Aft}) 1 * < Ai < •• • < Aat-i < 1 + (1 - A[t))-t 
ift is even and 

(2.9) 1 - (1 -A[t])t < Ai 
ift is odd. 

Theorem 3 (Bauer- Jost). Let B[t] be an upper bound for the largest eigenvalue 
Ajv_i[£] of A[t], i.e. Aat-i [t] < B[t}. (Since Aat-i [t] < 1 for t even, we can as- 
sume in this case w.l.o.g. that B[t] < 1 in this case.) Then all eigenvalues of A are 
contained in the union of the intervals 

o,i-(i-£[t])i]U[i + (i-£M)U 
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if t is even and 

Atf_i<l-(1-B[f])* 

if t is odd. 

These theorems show how random walks on graphs (or equivalently neighborhood 
graphs) can be used to estimate eigenvalues of the Laplace operator. In the rest of 
this paper we will use these insights to derive lower bounds for Ai and upper bounds 
for Ajv-i in terms of the Ollivier-Ricci curvature of a graph. 

2.3. Ollivier-Ricci curvature from a probabilistic view. We define a metric d 
on the set of vertices V as follows. For two distinct points x, y G V, d(x, y) is the 
number of edges in the shortest path connecting x and y. Then, including the family 
of probability measures m := {m x (-)}, we have a structure (V,d,m), on which the 
definition of Ricci curvature proposed by Ollivier [IS] can be stated. 

Definition 4 (Ollivier). For any two distinct points x,y G V, the (Ollivier-) Ricci 
curvature of (V, d, m) along (xy) is defined as 

Wx(m x ,m y ) 



(2.10) K(x,y):=l- 



d(x,y) 



Here, Wi(m x ,m y ) is the transportation distance between the two probability mea- 
sures m x and m y , in a formula, 

(2.11) W 1 (m x ,m y )= inf V V d(x', y')C' y (x', y'), 
where Yl{m x ,m y ) is the set of probability measures £ x,y that satisfy 

(2.12) t x ' v ( x '>y') = m *( x ')> E ^V.ifl = m M- 

y'ev x'ev 



The conditions (2.12) simply ensure that we start with the measure m x and end up 
with m y . Intuitively, Wi(m x ,m y ) is the minimal cost to transport the mass of m x to 
that of m y with the distance as the cost function. We also call such a £ x,y a transfer 
plan between m x and m y , or a coupling of two random walks governed by m x and m y 
respectively. Those (£ x,y might not be unique) which attain the infimum value 



in (2.11), are called optimal couplings. The optimal coupling exists in a very general 
setting. For locally finite graphs the existence follows from a simple and interesting 
argument in Remark 14.2 in [13] . 

There is a Kantorovich duality formula for transportation distances, 



(2.13) W x {m x ,m y )= sup 

/:Li P (/)<l 



fW)m x (x') - Y f(y')my(y') 
.x'ev y'ev 



where Lip(/) := sup^.^^ ■ For more details about this concept, we refer to 

Villani [221123], and Evans pf. 

Remark 3. Applying this definition of curvature to a Riemannian manifold with a 
family of probability measures obtained by restricting the volume measure to the closed 
r -balls, k will reduce to the Ricci curvature as r — >■ up to a scaling r 2 , see Ollivier 
[18] . for which we also refer for a general treatment of Ricci curvature. In Riemannian 
geometry, Ricci curvature controls how fast the geodesies emanating from the same 
point diverge on average, or equivalently how fast the volume of distance balls grows 
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as a function of the radius. Jost-Liu [TT] translate those ideas into a combinatorial 
setting and show that Ollivier-Ricci curvature on a locally finite graph reflects the 
relative abundance of triangles, which is captured by the local clustering coefficient 
introduced by Watts- Strogatz 



For the rest of this paper, let A; be a lower bound for the Ollivier-Ricci curvature, 



i.e. 



(2.14) K { x ,y) > k, Wx ~ y. 

Remark 4. By Proposition 19 in Ollivier [18j ; it follows that if k is a lower curvature 
bound for all neighbors x, y then it is a lower curvature bound for all pairs of vertices. 
This also follows from Theorem ^ below. 

Remark 5. By definition, the lower bound k for the curvature k is no larger than 
one. In fact, such a lower bound k always exists. Since the largest possible distance 
between points from the supports of two measures m x and m y at a pair of neighbors 
x,y is 3, we can easily estimate K,(x,y) > —2. We will derive a more precise lower 
bound for k on a locally finite graph with loops in Theorem^ see also Lin-Yau [1_5] 
and Jost-Liu [TT] for related results. 



We could also write (2.14) as 

(2.15) W 1 (m x , m y ) < (1 - k)d{x, y) = 1 - k, Vrr ~ y, 

which is essentially equivalent to the well known path coupling criterion on the state 
space of Markov chains used to study the mixing time of them (see Bubley-Dyer [3] or 
[13] . [T9"]). We will utilize this idea to interpret the lower bound of the Ollivier-Ricci 
curvature as a control on the expectation value of the distance between two coupled 
random walks. 

We reformulate Bubley-Dyer's theorem (see [3] or [13], [19]) in our language. 

Theorem 4 (Bubley-Dyer). On (V,d,m), if for each pair of neighbors x,y G V, we 
have the contraction 

Wi(m x , m y ) < (1 - k)d(x, y) = 1 - k, 
then for any two probability measures \i and v on V, we have 

W 1 (ixP,iyP)<(l-k)W l (ix,u). 

With this at hand, it is easy to see that if for any pair of neighbors x, y, k(x, y) > k, 
then for any time t and any two x, y, which are not necessarily neighbors, the following 
is true, 

(2.16) Wx&pt^F*) < (1 - kfd&y). 

We consider two coupled discrete time random walks (X t , Y t ), whose distributions are 
5 x P f , 5yP l respectively. They are coupled in a way that the probability 

p(X t = x\Y t = y')=^(x',y'), 

where Q' V {'-> ') i s the optimal coupling of S x P f and SyP 1 . In this language, we can in- 
terpret the term Wi(6 x P f , 5yP l ) as the expectation value of the distance Fj x '^d(X t , Y t ) 
between the coupled random walks X t and Y t . 

Corollary 1. On (V, d, m), if k(x, y) > k, Vx ~ y, then we have for any two x,y G V , 

(2.17) E^d{X t ,Y t ) = W^P'^yP 1 ) < (1 - kfd(x,y). 
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2.4. Ollivier-Ricci curvature and the lower bound of the first eigenvalue. 

In his paper [18] . Ollivier proved a spectral gap estimate which works on a general 
metric space with random walks. In particular, on finite graphs, it could be stated as 
follows. 

Theorem 5 (Ollivier). On (V,d,m), if n(x,y) > k, Vx ~ y, then the first eigenvalue 
Ai of the normalized graph Laplace operator A satisfies 

Ai > k. 

This is a discrete analogue of the estimate for the smallest nonzero eigenvalue 
of the Laplace-Beltrami operator on a Riemannian manifold by Lichnerowicz. As 
pointed out in Ollivier |18] . this result is also related to the coupling method for 
estimates of the first eigenvalue in the Riemannian setting developed by Chen- Wang 
[7] (which leads to a refinement of the eigenvalue estimate of Li-Yau [E]), see also 
the surveys Chen [SJ |6]. The corresponding result of Corollary [T] in the smooth case, 
i.e., controlling the expectation distance of two coupled Markov chains in terms of 
the lower bound of Ricci curvature on a Riemannian manifold, is a key step in Chen- 
Wang's method. 

A direct proof of Theorem [5] can be found in [TH] . Here we want to present an 
analogue of Chen- Wang's method in the discrete setting. Our proof of Theorem [5] 
motivated us to combine the Ollivier-Ricci curvature and the neighborhood graph 
method via random walks. It reflects the deep connection between eigenvalue esti- 
mates and random walks or heat equations. 

Proof: We consider the transition probability operator P : l 2 (V,fi) — > l 2 (V,fi) 
defined by Pf(x) := J2 y f(y) m x(y) = T, y fW x P(y)- Then we have P t f{x) = 
^2 f(y)S x P t (y). We construct a discrete time heat equation, 



(2.18) 



f(x,0) = h(x), 
f(x,l)-f(x,0) = Af(x,0), 
f(x,2)-f(x,l) = Af(x,l), 

f(x,t + l)-f(x,t) = Af(x,t), 



where fi(x) satisfies A/i(x) = — Ai/i(x) = Pfi(x) — fi(x). Iteratively, one can find 
the solution of the above system of equations as 

(2.19) f(x,t) = P t f 1 (x) = (l-X 1 ) t f 1 (x). 

We remark here that the solution of the heat equation on a Riemannian manifold with 
the first eigenfunction as the initial value is f(x,t) = fi(x)e~ Xlt , which also involves 
information about both the eigenvalue Ai and the eigenfunction fi(x). 

If Ai > 1 there is nothing to prove since k < 1 by definition. Therefore we can 
suppose Ai < 1 in the following. Then we have for any x,y G V 

(l-\Mi(x)-h(y)\ = \f{x,t)-f(y,t)\ 

= \pt fl (x)-p t f 1 {y)\ 

x' 

<Lip(/ 1 )E^(x t ,F t ) 
<Lip{f 1 ){l-k) t d{x,y). 
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Here, Lip(/) is always finite since the underlying space V is a finite set. In the last 
inequality we used Corollary [T} From an analytic point of view, the above calculation 
can be seen as a gradient estimate for the solution of the heat equation. 

Since the eigenfunction f\ for the eigenvalue Ai is orthonogal to the constant func- 
tion, i.e. (fx, 1) M = 0, we can always find x ,y E V such that \ fi(x ) — /i (2/0) I > 0. 
It follows that 



< ( — X J Lip(/i)d(x ,y )- 

To prevent a contradiction when t — > 00, we need 

1 - k 

(2.20) j-jj- > 1, 

which completes the proof. □ 

3. Estimates for Ollivier-Ricci curvature on locally finite graphs 

with loops 

Jost-Liu [H] give a sharp estimate for Ollivier-Ricci curvature on locally finite 
graphs without loops. As shown in Example [TJ neighborhood graphs are in general 
weighted graphs with loops. Therefore, for our purposes, we need to understand the 
curvature of graphs with loops. In this section, we generalize the estimates in Jost- 
Liu [TT] for locally finite graphs G = (V, E) that may have loops. This is done by 
considering a novel optimal transportation plan. 

We first fix some notations. For any two real numbers a, b, 

a + := max{a, 0}, a A b := min{a, b}, and a V b := max{a, b}. 

We denote N x := {z E V\z ~ x} as the neighborhood of x and N x := N x U {x}. 
Then N x = N x if x has a loop. For every pair of neighbors x, y, we divide A^, N y into 
disjoint parts as follows. 

(3.1) N x = {x} U {y} U N l x U N xy , N y = {y} U {x} U N] U N xy , 

where 

and 



N x ■= {A z ~ x,z ^ y,z ^ y}, 



TV T C I II XZ Z1] ~\ 

N x > y := {z\z ~ x,z ~ y,z ^ x,z TV,~f- > -p}, 

(l X Qjy 

TV T C I II ^-^XZ ZXI ~\ 

N x<y := {z\z ~ x, z ~ y, z ^ x, z ^ y, — < —}. 

Ct X Qjy 

In the next figure we illustrate this partition of the vertex set. 



1 

N x 




x<y 
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Theorem 6. On G = (V, E), we have for any pair of neighbors x, y e V, 
K(x,y)>k(x,y):=- [l-^-^- ^ ^ V W * lV 



y Xl£N X y 



J j ^sj/ ^oij/ W X1X ^ VJ Xll 

d x dy d x d y 



W XlX W Xl y W XX Wyy 

+ A A d d ' 

xl eN xy r y x y 

Moreover, this inequality is sharp. 

Remark 6. On an unweighted graph, the form k(x,y) for x ~ y becomes 

k( x y ) = - (i - — - — - M^yy\ _ (i _ JL _ JL _ M^yL\ 

' V d x d y d x Ad y J + \ d x d y d x V d y J + d x V d y d x d y 

where §(x,y) := Yl Xl eN xy 1 ^ s ^ e number of triangles containing x,y, c(x) = or 1 
is the number of loops at x. 



Proof : Since the total mass of m x is equal to one, we obtain from (3.1 ) the following 
identity for neighboring vertices x and y: 

/q q\ -, V^xy V ^ ^x±x 1U XX \ r t^xix 

[6 - Z) l ~^r~ 2^ ~a~ = ^r + 



dor. d x d x d x 

xi£N xy Xi6ATl 



A similar identity holds for y. 
We denote 



l^xy t^xy V v ^x\x W Xiy 



a _ i w xy w xy 
A x , y - 1 - -j- - -r- - 2^ 



V 



d x d y x ^ N ^ d x d y 

V^xy V^xy \ v 1^x\x . W xiy 



r> i _ w xy _ Wxy_ _ \ ^ 

'i d„ 2- d, 



"x wy 



xy 



Obviously, A x>y < B x<y . We firstly try to understand these two quantities. 
If A XjV > 0, we have 

^ ^ dy " ^ r dy ~ d,^ \ d X dy 

y x^n xv y x xi£N x > y v y 



i.e., using (3.2) we observe that the mass of m y at y and Ny is no smaller than that 
of m x at y and the excess mass at N x > y . Rewriting ( |3.3[ ) in the form 

W ^2/ _|_ ^ w xy / "W-'aa-i 3/ 

(l " x 1& N xy ^ dx r,( .V, ., V ^ 

and subtracting the term Y\ cW A on both sides we obtain 

o /-^ixi£Nxy d x d y 

UJ X n v v / VJ x ,'u UJxxi \ , - U)v 



X <y 



(3 4) xy _|_ / g ^ _ xx i j < x — _ 



i.e., the mass of m x at x and is larger than that of m y at x and the excess mass 



a t N x<y . 
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If B xy > 0, we have 



(3.5) 



w 



■>-y 



E 



XXl 



d, 



E 



w 



XXl 



Wxn 



> 



w 



i.e., the mass of m x at x and N% and the excess mass at N x > y is no smaller than that 
of m y at x. 

In Jost-Liu [11] it is explicitly described how much mass has to be moved from a 
vertex in N x to which point in N y , i.e. the exact value of £ x ' y (x', y'), for any x' G N x , 
y' G N y . But in the case with loops it would be too complicated if we try to do 
the same thing. Instead, we adopt here a dynamic strategy. That is, we think of a 
discrete time flow of mass. After one unit time, the mass flows forward for distance 
1 or stays there. We only need to determine the direction of the flow according to 
different cases. 

As in Jost-Liu [TTj . we divide the discussion into 3 cases. 

• < A x>y < B XjV . In this case we use the following transport plan: Suppose 
the initial time is t — 0. 

t — 1: Move all the mass at N* to x and the excess mass at N x > y to y. We 
denote the distribution of the mass after the first time step by m 1 . We 
have 



WAm^m 1 ) < 1 



W 



XI) 



d. 



d. 



xi£N xy 



d, 



x 1 + 



E 



w 



XXl 



W 



d, 



d„ 



x 1 



t = 2: Move one part of the excess mass at x now to fill the gap at N x<y 
and the other part to y. By (3.4) the mass at x after t = 1 is enough to 
do so. The distribution of the mass is now denoted by m 2 . We have 



£ E 

Xl&N X< y 



w 



W 



XXl 



d„ 



dx 



+ 



i _ 



E 

Xl&N X y 



X 1 



- E 

xi&N x<y 



W., 



xyy_ 

dy 



Wxy 

d„ 



x 1 



t = 3: Move the excess mass at y now to N*. We denote the mass after the 



third time step by m 



m y . We have 



Wi(m 2 ,m„) < 



E 

xi£N xy 



Wxxi 

dr 



- E 

%l£N X< y 



E 

Xl£N x >y 



w 



xxi 



d'n 



w 



yy 



X 1 



By triangle inequality and (2.11), we get 

Wx(m x , my) <Wi(m x , m l ) + W^m 1 , m 2 ) + W^m 2 , m 



1U Xl y 



d, 



=3-2 



w 



xy 



w 



xy 



d , 



d„ 



^ E 



w 



XX l 



d n 



A 



v) 

tU Xl y 
dy 



E 



w 



XXl 



w 



xy 



+ 



W 



d a 



d, 



V 



tUxn 

di, 



dn 



Xl£N X y y Xl£N X y 

Moreover, if the following function can be extended as a function on the graph 
such that Lip(/) < 1, (i.e., if there are no paths of length 1 between N* and 
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N x<y , nor paths of length 1 between Ny and N x > y , nor paths of length 1 or 2 
between and Ny,) 

if z G JVj; 

fi - , ■ • if z e M U iV x<y ; 
J{Z> \ o if z G {x} U N x > y ; 
if 2 G 



then by Kantorovich duality (2.13), we can show that the inequality above is 
actually an equality. Recalling the definition of n(x,y), we have proved the 
theorem in this case. 

A x>y < < B xy . We use the following transfer plan: 

t — 1: We divide the excess mass of at N x > y into two parts. One part 
together with the mass of m x at y is enough to fill gaps at y and Ny. 
Since (3.3) doesn't hold in this case, this is possible. We move this part 
of mass to y and the other part to x. We also move all the mass of m x 
at N* to x. 

t = 2: We move the excess mass at x now to N x<y and the excess mass at 
V to N}. 

Applying this transfer plan, we can prove (we omit the calculation here) 



W 1 (m x ,m v ) < 2 



w 



w 



xy 



d,, 



A 



■w. 



%iy 



d„ 



Wxx 

d x 



w 



d„ 



Moreover, if the following function can be extended as a function on the graph 
such that Lip(/) < 1, (i.e., if there are no paths of length 1 between N].UN x >y 
and N} U N x<y ,) 



0, if z G JVJ U N x<y - 

1, if z = x or z = y; 

2, if z G Nl U N x > y , 



then by Kantorovich duality (2.13), we can check that the inequality above is 
actually an equality. 

A X)V < B x>y < 0. We use the following transport plan: 



t — 1: Move the mass of m x at iV^ and N x > y to x. Since now (3.5) doesn't 
hold, we need to move one part of the mass m y (y) to x and the other 
part to Ny and N x<y . 
Applying this transfer plan, we can calculate 



W-\(m x ,m^ < 1 



E 



w 



xxi 



d, 



A 



w 



ui. 



w 



yu 



d, 



Since the following function can be extended as a function on the graph such 
that Lip(/) < 1, 



0, if z G {x} U N x<y U N^; 

1, if z G {y} U N x > y U N£, 



we can check the inequality above is in fact an equality by Kantorovich duality. 
That is, in this case for any x ~ y, 



K(x,y) 



E 



W; 



d 7 



A 



w 



xiy 



d„ 



+ 



Wxa 

d r 



+ 



w 



mi 



d„ 
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□ 

We also have a generalization of the upper bound in Jost-Liu [IT] on G. 
Theorem 7. On G = (V, E), we have for every pair of neighbors x, y, 

X!e{x}U{y}LlN xy X V 

Proof : / := Ylxie{x}u{y}uN x ~d^~ ^ ~^ JL * s exac tly the mass of m x which we need 
not move. The other mass need to be moved for at least distance 1. So we have 
Wi(m x , m y ) >1 — I, which implies k{x, y) < I, for x ~ y. □ 

Remark 7. If for every x 6 V, w xx = 0, then Theorem^ and Theorem^ reduce to 
the estimates in Jost-Liu 

Remark 8. Theorem^tells us that if k(x, y) > for any x ~ y, then either y) ^ 
or c(x) 7^ or c(y) ^ 0. 

Example 2. We consider a lazy random walk on an unweighted complete graph K.^ 
with N vertices governed by m x (y) = 1/N, Vx,y. Or equivalently , we consider the 
graph /C^ zy . Using Theorem^ and Theorem^ we get for any x, y 

N-2 11 . . 1 

That is, in this case, both the lower and the upper bound are sharp. 

4. Estimates of the spectrum in terms of Ollivier-Ricci curvature 

In this section, we use the neighborhood graph method to derive estimates for the 
spectrum of the normalized graph Laplace operator in terms of the Ollivier-Ricci 
curvature. In particular, we get a sharp upper bound for the largest eigenvalue using 
the curvature on the original graph. Furthermore, we give nontrivial estimates using 
the curvature on neighborhood graphs even if the curvature of the original graph is 
nonpositive. In this case, our results improve Theorem [5] 

Lemma 2. Let k[t] be a lower bound of Ollivier-Ricci curvature of the neighborhood 
graph G[t}. Then for all t > 1 the eigenvalues of A on G satisfy 

(4.1) 1 - (1 - k[t])* < Ai < • • • < Aat-i < 1 + (1 - k[t])* 
if t is even and 

(4.2) < A x 

if t is odd. Moreover, if G is not bipartite, then there exists a t' > 1 such that for all 
t > t' the eigenvalues of A on G satisfy 

< 1 - (1 - fc[t])* < Ai < • • ■ < Aiv-i < 1 + (1 - k[t])? < 2 
if t is even and 

< 1 - (1 - k[t])t < Ai 

if t is odd. 



Proof : By combining Theorem [2j and Theorem|5jwe directly obtain equation (4.1) 
and equation (4.2). 

The second part of this Lemma is proved in two steps. In the first step, we will show 
that if G is not bipartite then there exists a t' such that for all t > t' the neighborhood 
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graph G[t] of G satisfies w xy [t] ^ for all x, y G V, i.e. G[t] is a complete graph and 
each vertex has a loop. In the second step, we show that any graph that satisfies 
w xy 7^ for all x, y G V has a positive lower curvature bound, i.e. k > 0. This then 
completes the proof. 

Step 1: By the definition of the neighborhood graph it is sufficient to show that 
for all t > t' there exists a path of length t between any pair of vertices. Since G 
is not bipartite it follows that there exists a path of even and a path of odd length 
between any pair of vertices in the graph (This is the definition of a bipartite graph!). 
Given a path of length L between x and y then we can find a path of length L + 2 
between x and y as follows: We go in L steps from x to y and then from y to one of 
its neighbors and then back to y. This is a path of length L + 2 between x and y. 
Since G is finite, it follows that there exists a t' such that there are paths of length t 
for all t >t' for any pair of vertices. 

Step 2: Given a graph that satisfies w xy ^ for all x, y G V. 

Since each vertex in the graph is a neighbor of all other vertices, it is clear that we 
can move the excess mass of m x for distance 1 to anywhere. Therefore 

W 1 {m x ,m y ) < 1 - 2 , 

which implies 

K(x,y) > 



— A 
d x dy 



^ w xiy 



d x dy 



By Theorem [7j it follows that the above inequality is in fact an equality. Hence for 
all x, y G V, we have 

K ( x y\ _ ^ Wxx i /\ Wx ^V > ^y mm a;,y w xy > gim^ W xy ^ ^ 

^ e y ^ max x (ia, max XJ/ if^^ 

since the weight is positive for every pair of vertices. 

This completes the proof. □ 

Remark 9. From Remark^ after Theorem^ we know that positive lower curvature 
bound is a strong restriction on a graph. Hence, Olliver's estimate X± > k only yields 
nontrivial estimates for Ai in those restricted cases. Lemma shows that, unless G 
is bipartite, the neighborhood graph method can be used to always obtain a nontrivial 
lower bound for Ai in terms of the Ollivier-Ricci curvature. 

Now, we will show that it is also possible to control the largest eigenvalue of A on 
the graph G in terms of the Olliver-Ricci curvature of G itself. 

Theorem 8. On (V,d,m), if n(x,y) > k, Vx ~ y, then the largest eigenvalue Xn-i 
of the normalized graph Laplace operator A satisfies 

Ajv-i <2 — k. 

Remark 10. As discussed above, the largest eigenvalue Ajy-i is small if the graph 
does not contain a locally bipartite subgraph that has only a few connections to the 
rest of the graph. It is well known that a graph is bipartite if and only if it does 
not contain any cycle of odd length, in particular a bipartite graph does not contain 
triangles and loops. Since positive Ollivier-Ricci curvature implies that locally there 
always exists a triangle or a loop, i.e. the graph cannot be locally bipartite, it is not 
surprising that the largest eigenvalue can actually be controlled from above in terms 
of the Ollivier-Ricci curvature. 
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By combining Theorem [8] and Theorem |6j we have the following corollary. From 
Remark [6j we know that k(x, y) reflects the number of triangles and loops in the 
graph. 

Corollary 2. On G = (V,E), the largest eigenvalue satisfies 

Ajv-i < 2 — min k(x, y), 

where k(x,y) is defined in Theorem^ 

The main point in the proof of Theorem [8] is to explore the relation of the Ollivier- 
Ricci curvature k on the original graph G and the Ollivier-Ricci curvature K,[t] on its 
neighborhood graph G[t]. Before we prove Theorem [8] we consider some lemmata. 

If we interpret the graph G = (V,E) as a structure (V, d, m = {5 X P}), then 
by ( 2.7[ ) its neighborhood graph G[t] = (V, E[t]) can be considered as a structure 
(V, d[t], {8 X P 1 }). So the first step should be to estimate the variance of the metrics 
on neighborhood graphs. 

Lemma 3. For any x, y G V, we have 

(4.3) ^d{x,y)<d[t]{x,y). 

Proof : For any x, y G V, we set c?[t](x, y) = oo if we cannot find a path connecting 
them in G[t\. Otherwise, we just choose a shortest path xq = x, xi, . . . , x\ = y, 
between x and y in G[t], i.e. / = <i[t](x, y). For Xi, Xj+i, i = 0, — 1, by definition 
of neighborhood graph, we have d(xi,Xi + i) < t in G. Equivalently, 

^d(xi,x i+1 ) < 1 = d[t](xi,x i+1 ). 
Summing over all i, we get 

2-1 



^-^2d(xi,x i+1 ) <d[t](x,y). 



Then the triangle inequality of d on G gives (4.3). □ 

Remark 11. In fact, when t is larger than the diameter D of the graph G, we have 
a better estimate 

(4.4) \d{x,y) < ^d(x,y) < 1 < d[t]{x,y). 

Lemma 4. If E C E[t], then d[t](x,y) < d(x,y). 

The proof is obvious. The interesting point is that when the Ollivier-Ricci curvature 
of the graph G is positive, E C E[t] is satisfied for all t and hence Lemma [1] is 
applicable. This can be seen as follows. First we observe that for (x, y) G E, if 
§(x,y) 7^ or c(x) ^ or c(y) ^ 0, then (x,y) G E[t]. Then we get the conclusion 
immediately from Remark [8] after Theorem [7j 

Lemma 5. Let k be a lower bound of k on G. If E C E[t], then the curvature n[t] 
of the neighborhood graph G[t] satisfies 

(4.5) K[t](x,y) > l-t(l-Jfe)*, \/x,yeV. 
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Proof : By Lemma |4j Corollary [T] and Lemma |3j we get 

wf* ] (<y^*, V*) < <(<U>*AP') 

< (l-kYd(x,y) 

< t(l - kYd[t}(x,y). 

We use wf^, Wf here to indicate the different cost functions used in these two 
quantities. In the first inequality above we used that the transportation distance 
(2.11) is linear in the graph distance d(-, •). Recalling the definition of the curvature, 
we have proved (4.5). □ 
Now we arrive at the point to prove the upper bound of the largest eigenvalue. 
Proof of Theorem [§} Using Lemma [5] and Theorem [5j we know on G[t), 

Ai[t] > l-*(l-jfe)*. 

Then by Lemma [2j we get for any even number t, 

Xn-i < 1 + ^(1 — k). 

Letting t — > +oo, we get \n-i < 2 — k. □ 
By combining Theorem [5] and Theorem [8j we get 

k < Ai < • • • < Ajv-i <2 — k, 



That is, (4.1) is also true for t = 1. 

Example 3. On an unweighted complete graph JC^ with N vertices, we have 



JV-l' 1 Af-1 

Therefore, 

k < Ai = • • • = Aat_i = 2 — k. 
That is, our upper bound estimate for is sharp for unweighted complete graphs. 

Example 4. Let's revisit the graph /C^ zy in Example^ We have 
k = k(x, y) = 1, Vx, y and A = • • • = Ajv_i = 1. 

Therefore, 

k — Ai = • • • = Xn-i = 2 — k = 1. 
That is, both estimates are sharp in this case. 



In fact, we can show that (4.1 ) is true for any t > 1. 

Theorem 9. Let k[t] be a lower bound of Ollivier-Ricci curvature of the neighborhood 
graph G[t]. Then for all t > 1 the eigenvalues of A on G satisfy 

(4.6) 1 - (1 - fc[*])i < Ai < • • • < \ N -! < 1 + (1 - fc[t])i 

Moreover, if G is not bipartite, then there exists a t! > 1 such that for all t > t' the 
eigenvalues of A on G satisfy 

< 1 - (1 - k[t]) l t < Ai < • • • < \ N -! < 1 + (1 - k[t])? < 2. 

Proof : By LemmajiJ we only need to prove the upper estimate Aat_i < 
if t is odd. This follows immediately from Theorem [8] and Theorem |3j □ 

In particular, unless G is bipartite, we can obtain nontrivial lower bounds for Ai 
and nontrivial upper bounds for \n~i even if k = k[l] is nonpositive. 
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5. Estimates for the largest eigenvalue in terms of the number of 

joint neighbors 

In Bauer- Jost [1] it is shown that the next lemma is a simple consequence of The- 
orem [TJ 

Lemma 6. Let u be an eigenfunction of A for the eigenvalue A. Then, 

(u,A[2]u), _ Y, x , y w X y[2}{u{x)-u{y)f 



(5.1) 2 - A 



(w,Aw) M J2 x , y Wxy(u(x)-u(y)) 2 



Lemma [6] can be used to derive further estimates for the largest eigenvalue Aat__i 
form above and below. We introduce the following notations: 

Definition 5. Let N x be the neighborhood of vertex x as in Section [i| The mini- 
mal and the maximal number of joint neighbors of any two neighboring vertices is 
defined as | a := min^f^x, y) + c(x) + c(y)) and \ 2 : = max^f^x, y) + c(x) + c(y)), 
respectively. Furthermore, we define W := max Xiy w xy and w := m.in Xjy]X ^ y w xy . 

Theorem 10. We have the following estimates for Ajv-i- 
(z) If E{G) C E{G[2]) then 

Ajv-i < 2 



W maXa; d x 
(it) IfE(G[2]) C E[G) then 

„ W 2 h ^ x 

S Aat_i 



w min r oL 



Proof. On the one hand, we observe that if E(G) C i£(G[2]) then for every pair of 
neighboring vertices x ~ y in G 

(52) ^ = ?i^ > 



w 2 



w I9 w xy W max x d x 

On the other hand if E(G[2]) C -E(G) then for every pair of neighboring vertices 
x ~ y in G(2) we have 

(5 3) Wxy [2] = E z T z W xz W zy < jj 2 



u?^ it; min^. 



Substituting the inequalities (5.2) and (5.3) in equation (5.1 ) completes the proof. □ 



Interestingly, Theorem 10 (z) yields an alternative proof of Theorem |8| in the case 
of unweighted regular graphs. Since Theorem [8] trivially holds if k < we only have 
to consider the case when k > is a lower curvature bound. The discussion after 
Lemma |4| sh ows that k > implies that E(G) C J5(G[2]) and hence we can apply 



Theorem 10 (i) in this case. From Theorem \f\ it follows that for an unweighted graph 



<*,y)<^Mr + c{x) ' c(y) 



\/ ditf ^ Jr ti 



for all pairs or neighboring vertices x, y. In the case of an <i-regular graph G this 
implies that a lower bound k for the Ollivier-Ricci curvature must satisfy, 
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Hence for an unweighted ci-regular graph Theorem 10 implies 



Ajv-i <2-% <2-k. 
a 

We consider the following example. 

Example 5. Let's revisit the complete unweighted graph JC^. For the complete graph 
we have E(JCn) C E(JCn[2]). We have fji = N — 2 and max x d x = N — 1. Thus, 
Theorem^T^(i) yields 

\ N 
Xn -' ~ iV^T' 

i. e. the estimate from above is sharp for complete graphs. Now we consider the 
unweighted complete graph on N vertices with N loops /C^ zy . We have E(K ^ y \2] ) 
E(K}^ Y ). Furthermore, we have $2 = N and mm x d x = N. Thus, Theorem 
yields 

i. e. the estimate from below is sharp for complete graphs with loops. 
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6. An example 

In this section, we explore a particular example, the circle C5 with 5 vertices. We 
show that our estimates using the neighborhood graph method can yield nontrivial 
estimates although the curvature of the original graph has a nonpositive lower curva- 
ture bound. We also discuss the growth rate of the lower bound k for curvature K[t] 
as t — > 00 on C5. 

We consider the unweighted graph C 5 displayed in the next figure. 




We know that the first and largest eigenvalue of A on C5 are given by 

Ai = 1 - cos — = 0.6910, A 4 = 1 - cos — = 1.8090. 
5 5 

It is easy to check that the optimal lower bound k for the curvature is 0. So in this 
case Ollivier's result Theorem [5] and Theorem [8] only yield trivial estimates. 

Now we consider the neighborhood graph C§[2] depicted in the next figure (we 
change the order of vertices). 
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The weight of every dashed loop is 1 and the weight of every solid edge is 1/2. We can 
check that the optimal k[2] is 1/4. Then Theorem M yields the nontrivial estimate, 



AT /o 

Ai > 1 - — = 0.1340, A 4 < 1 + — = 1.8660. 
2 ' 2 



Moreover, the neighborhood graph Cs[4] = (C[2])[2] is depicted in the next figure. 




The weight of every dashed loop is 3/4 and the weight of every solid edge is 1/2 and 
every dash-dotted edge is 1/8. We can check that the optimal lower curvature bound 
is k [4] = 1/2. So Theorem M tells us that 



Ai > 1 - 4t= = 0-1591, A 4 < 1 + 4= = 1-8409. 
</2 #2 



We can also consider the neighborhood graph of odd order Cs[3], which is depicted 
in the next figure. 
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The weight of every solid edge is 3/4 and the weight of every dashed edge is 1/4. 
Then the optimal lower curvature bound is given by fc[3] = 3/8. Theorem M implies 



Ai > 1 



= 0.1450, A 4 < 1 



1.8550. 



From the above calculations, we see that the Ollivier-Ricci curvature n[t] on the 
neighborhood graphs give better and better estimates in this example. 



We also further compare the estimates in Theorem 10 and Theorem M on the 



weighted graph C 5 [3]. By Theorem 10 we get 



A 4 [3] < 2 - 



(V4) a 
3/4 



3 
2 



15 



13 

S 



This is worse than the one given by Theorem 8l A 4 [3] < 2 — j 

Another interesting problem is concerned with the limit of the neighborhood graphs. 
As shown in Bauer- Jost pQ, since C5 is a regular non-bipartite graph, C 5 [t] will converge 
to C 5 [oo] := C 5 as t — > 00. In this case C 5 is a complete graph and every vertex has 
a loop and the weight of every edge in C 5 is 2/5. We can then check that k[oo] = 1. 
By Theorem [91 for large enough t, we have 



(6.i; 



< 1 - (1 - k[t))t < Ai 



27r 1 
1 — cos — < 1. 

5 



Therefore < (1 — k[t])* < 1. Intuitively, k[t] should become larger and larger as 
t — > 00 since the graph C 5 [t] converges to a complete graph and its weights become 
more and more the same. We suppose lim t _ ) . 0O (l — k[t])t exists (this is true at least 
for a subsequence). Then to avoid contradictions in ( |6.1[ ), we know there exists a 
positive number a such that 



lim(l - k[t]Y 

t— >co 



That is, 



lim 

t— >oo 



log(i-fcM) 



> 0. 



—a. 



which means k[t] behaves like 1 — P(t)e at as t — > 00 where P(t) is a polynomial in t. 
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